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This study is designed to explore the two dimensional temperature maps of the atmospheric argon discharge 
consisting of pin-to-plane electrodes supplied by a high voltage DC source. After checking the stability of 
the micro discharge, the two dimensional image plane focused by a quartz lens was scanned by the fiber 
probe driven by a 3D Mobile Platform. The rotational and vibrational temperatures are calculated using 
nitrogen emissions collected by the high resolution spectrometer and high sensitive intensified charge coupled 
device (ICCD). The rotational temperature varies from 1558.15 K to 2621.14 K and vibrational temperature 
varies from 3010.38 K to 3774.69 K, indicating a great temperature gradient due to small discharge size. The 
temperature maps show a lateral expansion and a sharp truncation in the radial direction. A double layers 
discharge is identified, where an arc discharge coats the glow discharge. 

PACS numbers: 52.70.-m, 52.30.-q, 52.80.Hc, 52.80.Mg 


I. INTRODUCTION 

Atmospheric pressure plasmas have received much at¬ 
tention as a promising source in numerous applications 
recentlyThey can be used as anatytical tools, pho¬ 
tonic devices, in the surface treatment of vulnerable ma¬ 
terials, bacterial inactivation processes, surface modifi¬ 
cation, chemical vapor deposition, plasma etching on 
films^”^^. 

Characteristics of micro plasmas, such as neutral gas 
temperature, electron temperature, electron density, ion 
density, electric field and so on, play an important role 
in physical and chemical processes^^”^^. There are many 
studies concerning the relationship between the bulk gas 
temperature and control parameters, for example, dis¬ 
charge current and gas flux, where the gas tempera¬ 
ture is estimated using radiative transitions of diatomic 
molecules^^”^^. However, the temperature mentioned 
above is the mean temperature without spatial resolu¬ 
tion. Papers seldom reported the spatial resolution tem¬ 
perature maps of micro plasmas (micro discharges) al¬ 
though which can give more information of micro dis¬ 
charge mechanism and help optimization of micro dis¬ 
charge device for applications. 

In this work, an experiment was designed to obtain the 
two dimensional gas temperature maps of micro plasmas 
generated in the atmosphere of the ambient air between 
a miniature flow argon and a metal film. Firstly voltage- 
current curve was measured. Then the stability of the 
micro discharge was checked and image plane focused by 
a quartz lens was scanned by a fiber probe and two di¬ 
mensional spectrums of nitrogen emissions were got by 
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a high resolution spectrometer and an intensified charge 
coupled device (ICCD). At last we got the two dimen¬ 
sional temperature maps of micro discharges suggesting 
a double layers coated discharge. 


II. EXPERIMENT 

The schematic of the experimental device is shown in 
Fig. 1. Argon gas flow jets into the ambient air through 
a stainless-steel capillary acting as the anode, whose in¬ 
ternal diameter is 175 iim and external diameter is one 
sixteenth of an inch. A metal film with smooth rounded 
edges serves as the cathode. The gap distance from the 
exit of capillary to the metal film is 2 mm. The circuit 
is driven by a high-voltage DC source with a ballast re¬ 
sistor limiting the discharge current. The argon gas flux 
is set as 10 seem controlled by a mass flow controller. So 
the mean speed of gas jet at the exit nearly equals 6.93 
m/s. The titanium film used in this work is achieved by 
low temperature deposition on silicon substrate using re¬ 
active DC magnetron sputtering. More details about the 
deposition equipments can be found in the referance^^. 

The primary methodology to get two dimensional tem¬ 
perature maps is that fiber probe driven by a 3D Mobile 
Platform scans the image plane focused by a quartz lens. 
Because signals of spontaneous radiation collected by the 
fiber probe is very weak. A high sensitive and enhance¬ 
ment ICCD (intensified charge coupled device, PI-MAX 
4) is used to read the spectrum data made by a high 
resolution spectrometer (2400 groove/mm grating). 

The typical DC discharge image is taken by a Manta 
G201C CCD with commercial image lenses at the dis¬ 
charge current of 10 mA shown in Fig. 2 (a). The visi¬ 
ble discharge structure is similar to that of a typical DC 
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FIG. 1. Schematic diagram of the experimental device, not 
scaled. 



FIG. 2. (a) Discharge image taken by the Manta G201G GGD 
with the commercial lenses, (b) The volt age-current curves of 
micro discharges. The error bars are the standard deviation 
of six measurements, (c) Air density distribution obtained 
by solving full component Naiver-Stokes equations, (d) the 
sketch map of radiation projection. 


glow discharge, namely a negative glow, a Faraday dark 
space and a positive column are observed. Furthermore, 
Fig. 2 (b) gives the voltage-current curve, which is the 
average of six measurements. The error bars are the stan¬ 
dard deviations of six measurements. The discharge volt¬ 
age drops slowly as the discharge increases suggesting a 
characteristic of arc discharge different from the glow dis¬ 
charge. The discharge image and voltage-current curve 
indicate the complex nature of the argon discharge. 

In this study, the second positive systems of nitrogen 
were used to estimate the rotational and vibrational tem¬ 
peratures assuming Boltzmann distribution of rotational 


and vibrational states population. We had better learn 
knowledge of the flow-field and air spatial distribution be¬ 
fore discharge. The full component Naiver-Stokes equa¬ 
tions were solved to estimate the air or nitrogen density 
distribution. The governing conservation equations are 
written below in cartesian coordinates: 


dU dE dF 
dt ^ dx ^ dz 
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and p is the mixture gas density, v are th e x and 2 : 
components of velocity respectively and V = is 

the speed, Et = e + ^ is the mixture gas energy density, 
Txx and Tzz are the mixture gas normal viscous stresses, 
Txz and Tzx are the mixture gas shear stresses. Ah = 
{cpi — Cp 2 )T is the enthalpy difference between argon and 
air, Cpi and Cp 2 are specific heat at constant pressure 
of argon and air respectively. D is the binary diffusion 
coefficient, ci and C 2 are argon and air mass fraction, 
C 2 = 1 — Cl. fx and fz are the x and 2 : components of 
volume force. 

Fig. 2 (c) shows the air density distribution. The air 
occupied the whole space until turning on the argon sup¬ 
ply. The argon gas spread outwards from the exit and 
brought an argon channel while the air was put away. 
However, the argon took over the argon channel rather 
than the whole space. The more further away the ar¬ 
gon channel, the more air or nitrogen. Fig. 2 (c) just 
shows the prominent aspect of the air distribution in two 
dimensional. In fact, jet flow and discharge are three di¬ 
mensional. Even x-z plane spatially resolved, the light 
on the image plane is the projective signals in the y di¬ 
rection. The projection geometry is shown in Fig. 2 (d). 

After knowing the projective nature of radiation, the 
spatial resolution of two dimensional spectrum depends 
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FIG. 3. Saturated exposure discharge images, (a) discharge 
image taken by the Manta G201C GGD with the quartz lens, 
(b) discharge image taken by the Manta G201G GGD with 
the commercial lenses. 


on the quality of picture on image plane focused by the 
quartz lens. Fig. 3 (a) shows the discharge image on im¬ 
age plane directly taken by the Manta G201C CCD with 
a high gain number and long exposure time, that is to 
say, a saturated exposure image. Different from Fig. 2 
(a), an obvious halo was found and the faraday dark 
space became unclear. At same time, a saturated ex¬ 
posure image was also taken by the Manta G201C CCD 
with the commercial lenses shown by the Fig. 3 (b). A 
coated discharge was identified where the discharge coat 
(coat discharge) surrounds the discharge body (body dis¬ 
charge). The nature of the coat discharge and the body 
discharge will be discussed in the future. 

Because two dimensional spectrums were collected 
non-simultaneously, the stability of spectrum should be 
checked. The continuous spectrum from 392 nm to 400 
nm was considered, which is emitted by = 

2 ^ = 5) and N 2 {C^Ilu.v' = 1 ^ B^Tlg.v" = 

4) where B^Ilg denote electron configuration, 

v'^ v" denote vibrational states. Emission intensity of 
rotational and vibrational bands is proportional to spon¬ 
taneous emission coefficient and population of up state, 

- AyfJlyllJN^ ( 3 ) 

where Ayfjfynjff is spontaneous emission coefficient, 

Ny>j> is population of up state, Xv'J'v"J" is wavelength, 
v' J' denote the upper vibrational and rotational state, 
v" J" denote the lower vibrational and rotational state, h 
is the planck constant, c is the speed of light. Under the 
Born-Oppenheimer approximation and Boltzmann distri¬ 
bution assumption, the intensity can be written as 

IbAj" = (4) 

where D is the constant, qv>,v" is the Franck-Condon fac¬ 

tor. Sj'^j" is the line strength intensity. Eyf and Ej/ are 
vibrational state energy and rotational state energy. The 
more details can be found in references^^’^^’^^’^^. 

Fig. 4 (a) shows curves of the rotational temperature, 
vibrational temperature and the intensity versus time of 
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FIG. 4. (a) curves of rotational, vibrational temperature and 
counts number versus time, (b) rotational and vibrational 
temperature distribution in the z direction, (c) lateral rota¬ 
tional and vibrational temperature distributions at z = 0. 


the test point on image plane. The intensity is the ra¬ 
diation intensity from 392 nm to 400 nm and recovered 
from the spectral data: 


Intensity{y^ z) 


Counts{y, z) 
time X gain 


( 5 ) 


where time is the integration time and gain is the amplify 
gain. 

The averages of rotational and vibrational tempera¬ 
tures are 2003.08 and 2770.32 K respectively. And the 
ratios between the standard deviation to the average of 
rotational temperature and vibrational temperature are 
4.60% and 1.64% respectively. The rotational tempera¬ 
ture and vibrational temperature do not depend on time, 
so the micro discharge can be regarded as stable. 

Fig. 4 (b) shows the one dimensional temperature dis¬ 
tribution along the axis of the steel capillary in the z 
direction. The rotational temperature and vibrational 
temperature decrease in the Faraday dark space, grow up 
in the positive column region and then decreases again. 
Fig. 4 (c) shows the lateral one dimensional temperature 
distribution at the z = t). It is found that the outer 
temperature is higher than inter temperature. 
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Fig. 5 (a) and Fig. 5 (b) show the two dimensional dis¬ 
tribution of rotational temperature and vibrational tem¬ 
perature, respectively. The rotational temperature varies 
from 1558.15 K to 2621.14 K and vibrational temperature 
varies from 3010.38 K to 3774.69 K. So a great temper¬ 
ature gradient exits in the small size discharge, which 
maybe produce large heat exchange. In addition. Fig. 5 
(c) shows the ratio map between rotational and vibra¬ 
tional temperature suggesting the disequilibrium degree 
of discharge. Fig. 5 (d) shows the discharge image on 
wave band from 392 nm to 400 nm recovered from spec¬ 
tral data in square root scale. The most significant fea¬ 
tures of temperature maps are the sharp truncation and 
the lateral expansion. The sharp truncation of temper¬ 
ature maps indicate a sudden interface between micro 
discharge and ambient air. In other words, there is not a 
continuous transition region connecting the ambient air 
and the micro discharge. 

The Fig. 5 (d) shows that intensity decreases rapidly 
in the radial direction, but the signals are still detected 
by the ICCD in the outer region resulting from the coat 
discharge. So the lateral expansions are observed in tem¬ 
perature maps. The outer temperature is higher than 
inter temperature, as noticed, the voltage-current curve 
suggested an arc discharge. Therefore, we can infer that 
the body discharge is the glow discharge while the coat 
discharge is arc discharge. 

The air or nitrogen density distribution shown in Fig. 2 
(c) brings a question why nitrogen emission image looks 
like Fig. 5 (d). According to Eq. 4, under the measured 
temperature distribution, the outer intensity should be 
stronger than inner intensity (obviously shown by Fig. 3 
(c)), as opposed to Fig. 5 (d). How does nitrogen disap¬ 
pear or look like less. A more possible reason is that little 
nitrogen molecules participate in the coat discharge. The 
coat discharge may be in filamentary mode unlike the dis¬ 
persive body discharge. Further more, the main excita¬ 
tion sources of nitrogen are electron, excited argon atom 
and ion. Besides electron excitation process, the heat 
transport through collision between nitrogen molecules 
and excited argon atoms may play an more important 
role. As we see from the temperature maps, the temper¬ 
ature is not maximum at the same x location, but the 
argon density and intensity are maximum. The energy 
transport of plasma particles is very complex, so the ef¬ 
fect of collision excitation processes should be explored 
by detailed simulation. 


III. CONCLUSION 

In conclusion, experiments were designed to get the 
two dimensional spectrums of the micro atmospheric ar¬ 
gon discharge on image plane collected by the high resolu¬ 
tion spectrometer and high sensitive ICCD. And then the 
two dimensional temperature maps of the micro plasma 
are obtained. The rotational temperature varies from 
1558.15 K to 2621.14 K and vibrational temperature 




FIG. 5. (a) Two dimensional map of the rotational temper¬ 
ature. (b) Two dimensional map of the vibrational tempera¬ 
ture. (c) Two dimensional map of ratio between the rotational 
temperature and the vibrational temperature, (d) The square 
root of discharge intensity on the wave band from 392 to 400 
nm recovered from spectral data. 


varies from 3010.38 K to 3774.69 K. So a great tem¬ 
perature gradient exits in the small size discharge, which 
will produce large heat exchange. The temperature maps 
show a lateral expansion and a sharp truncation in the 
radial direction. The temperature maps and saturated 
exposure discharge images suggest that a filamentary arc 
coat discharge surround the body glow discharge. And 
the double layers coated discharge is identified. 

The saturated exposure discharge images show a clear 
edge between the body and coat discharge, which is not 
reflected by the temperature maps. More higher spatial 
resolution researches are waited to study. And the two di¬ 
mensional electron temperature maps is being prepared. 


ACKNOWLEDGMENTS 


XXZ and SFX acknowledge support from the NSFC 
(Grant No. 11275127, 90923005), STCSM (Grant No. 
09ZR1414600) and MOST of China. 

^S.-Z. Li, J.-P. Lim, J. G. Kang, and H. S. Uhm, Physics of 
Plasmas (1994-present) 13 , 093503 (2006). 

^H. Shi, Y. Wang, and D. Wang, Physics of Plasmas (1994- 
present) 15 , 122306 (2008). 

^S.-Z. Li, Q. Wu, W. Yan, D. Wang, and H. S. Uhm, Physics of 
Plasmas (1994-present) 18 , 103502 (2011). 

^K. Ostrikov, Rev. Mod. Phys. 77 , 489-511 (2005). 

^K. K. Ostrikov, U. Cvelbar, and A. B. Murphy, Journal of 
Physics D: Applied Physics 44, 174001 (2011). 

®K. Ostrikov, E. C. Neyts, and M. Meyyappan, Advances in 
Physics 62 , 113-224 (2013). 











5 


^M. Biscan, Z. Kregar, U. Cvelbar, M. Mozetic, S. Milosevic, 
et al, Spectrochimica Acta Part B: Atomic Spectroscopy 103 , 
124-130 (2015). 

®S. J. Pearton, C. R. Abernathy, and F. Ren, Applied Physics 
Letters 64 (1994). 

®R. M. Sankaran and K. P. Giapis, Applied Physics Letters 79 

( 2001 ). 

^®P. B. Lukins, M. H. Zareie, and J. Khachan, Applied Physics 
Letters 78 (2001). 

^^B. Wu and A. Kumar, Applied Physics Letters 90 , 063105 (2007). 

^^P. J.-B. L. J.-H. L. J.-T. Kyung, Se-Jin, Applied Physics Letters 
91 , 091504 - 091504-3 (2007). 

Deng, C. Leys, D. Vujosevic, V. Vuksanovic, U. Cvelbar, 
N. De Geyter, R. Morent, and A. Nikiforov, Plasma Processes 
and Polymers 11 , 921-930 (2014). 

^^G. M. Petrov and D. Zhechev, Physics of Plasmas (1994-present) 
9 (2002). 

^^B. Qi, J. Huang, L. Gao, and Y. Qiu, Physics of Plasmas (1994- 
present) 16 , 083301 (2009). 


^®B. Qi, L. Pan, Q. Zhou, J. Huang, and Y. Liu, Physics of Plasmas 
(1994-present) 21 , 123302 (2014). 

^^G. Faure and S. M. Shkol’nik, Journal of Physics D: Applied 
Physics 31 , 1212 (1998). 

^®B. A. Cruden, M. V. V. S. Rao, S. P. Sharma, and M. Meyyap- 
pan. Applied Physics Letters 81 (2002). 

Staack, B. Farouk, A. Gutsol, and A. Fridman, Plasma 
Sources Science and Technology 14, 700 (2005). 

^^B. Bai, H. H. Sawin, and B. A. Cruden, Journal of Applied 
Physics 99, 013308 (2006). 

21y. Lu, S. F. Xu, X. X. Zhong, K. Ostrikov, U. Cvelbar, and 
D. Mariotti, FPL (Europhysics Letters) 102, 15002 (2013). 

^^S. F. Xu and X. X. Zhong, Applied Physics Letters 103 , 024101 
(2013). 

23w. Zhou, X. Zhong, X. Wu, L. Yuan, Q. Shu, W. Li, and Y. Xia, 
Journal of Physics D: Applied Physics 40, 219 (2007). 

M. Phillips, Journal of Physics D: Applied Physics 9, 507 
(1976). 

^^G. Hartmann and P. C. Johnson, Journal of Physics B: Atomic 
and Molecular Physics 11 , 1597 (1978). 



